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Vortex state structure of a Bose condensate in an asymmetric trap 
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Abstract 

Based on an analytic solution of the Gross-Pitaevskii equation in the large-condensate (Thomas-Fermi) limit 
we determine the structure of a stationary vortex in a Bose-Einstein condensate in a nonaxisymmetric rotating 
trap. The condensate velocity field has cylindrical symmetry only near the vortex core and becomes intrinsically 
anisotropic near the condensate boundaries. Rotating the anisotropic trap induces an additional irrotational 
velocity field even for a vortex-free condensate. 
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The experimental achievement of Bose-Einstein 
condensation in confined alkali-atom gases [1-3] 
has stimulated great interest in the generation and 
observation of vortices in such systems [4-6] . Ro- 
tating a totally anisotropic harmonic trap at an an- 
gular frequency f2 can, in principle, generate vor- 
tices; they are energetically stable for Q > £l c [7- 
10]. ^ 

This work determines the condensate veloc- 
ity field for a stationary vortex when the com- 
pletely anisotropic harmonic trap potential Vt r = 
\M{uj 2 x 2 +ui 2 y 2 +u 2 z 2 ) rotates about a principal 
axis. We consider a zero-temperature condensate 
in the Thomas-Fermi (TF) limit, when the vortex 
core radius £ ~ d 2 /R is small compared to the 
mean oscillator length d and the mean dimen- 
sion R of the condensate [here, d = {d x d y d z ) l ' s 



with di = y/h/Midi and trap frequencies u>i (i — 
x,y,z)]. 

We assume that the trap rotates with an angu- 
lar velocity Q around the z axis. At zero tempera- 
ture in the rotating frame, the equilibrium conden- 
sate wave function 'F satisfies the Gross-Pitaevskii 
equation: 

t-2 
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—V 2 + Vtr + gW - ix + imd^j * = 0,(1) 

where g — Airh 2 a/M > is the effective interparti- 
cle interaction strength, ft is the chemical potential 
in the rotating frame, and cf> is the azimuthal angle 
in cylindrical polar coordinates (r, cf>, z). Consider 
a g-fold vortex, oriented along the z axis, with a 
condensate wave function 'F = |\&|e . Its phase 
must have the form S = qcj> + Si , where Si is an 
odd 27r-periodic function of (/>. 

Multiplying the imaginary part of (1) by |^/| 
yields the equation of current continuity in the ro- 
tating frame: 
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div(|*| 2 VS) -^Sld |*| 2 = O. (2) 

The second term here arises from the difference 
between the condensate velocity field v' = jj^S 
in the laboratory frame and v r = -A VS — 51 x r in 
the rotating frame. In the TF limit in the absence 
of a vortex, g\^TF\ 2 = M — Hn and for distances far 
from the vortex core (r 3> \q\Ct we may substitute 
\^tf\ into (2) instead of |*|, getting 

^TF\ 2 V 2 S-[LU 2 x + Lul + (u;l-u; 2 y )co S (2^}rdrS 

^r 2 )-2u J 2 zd z S = 0. 



+(^-^)sin(20)l^5 

Let us seek a solution of this equation in the 
form: 



S = S + qcj) - /30r 2 sin(20) 



(3) 



where (3 = (M/2K) (w 2 - w 2 ) / (cj 2 + lo 2 ) . The re- 
sulting equation for So does not contain Q. Fur- 
ther, let us consider distances x,y,z <C R x , R V} R Z , 
where Ri are the TF dimensions of the condensate 
{MuSfR 2 = 2/j, i = x,y,z). At these intermediate 
distances we can omit terms of higher order in the 
parameter (xi/Ri) . As a result, we find the fol- 
lowing inhomogeneous equation: 

V 2 So+(^-^)«sin(20)=O . (4) 

This equation has an explicit periodic solution 



So = -aqr 2 In ( — j sin(20) , 



(5) 



where a — \ (R x 2 — Ry 2 ) an d ^4. is a constant of 
integration that must be chosen to minimize the 
total energy. We estimate A ~ Rx, where R\ = 
2i? 2 i? 2 /( J R 2 +R 2 y ). Thus, the phase S near the 
center of a rotating trap in the presence of a vortex 
is 



S: 



aqhi 



R~ ] + ' <l 



r 2 sin(2</>) . (6) 



Near the vortex core, the condensate wave func- 
tion and condensate velocity possess cylindrical 
symmetry, while far from the vortex core the con- 
densate velocity adjusts to the anisotropy of the 



trap and becomes asymmetric. The condensate ve- 
locity in the laboratory frame becomes 
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rcos(2</>) > e^ 



(7) 



The energy barrier to vortex formation [9] means 
that the condensate can be vortex free even for 
O > £7 C , in which case the trap rotation induces 
an irrotational velocity field v' oc il\7(xy) [11]. 
To logarithmic accuracy, the additional phase So 
does not affect the critical angular velocity £l c and 
metastable angular velocity Vl m = |O c found pre- 
viously [9]. 

In conclusion, an analytic solution of Gross- 
Pitaevskii equation in the rotating frame gives 
the condensate velocity field for a vortex state in 
an anisotropic rotating trap. Even for a vortex- 
free condensate, the rotation induces superfluid 
motion. 
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